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Abstract Issues concerning tlic structure and evolu- 
tion of core collapse progenitor stars are discussed with 
an emphasis on interior evolution. We describe a pro- 
gram designed to investigate the transport and mixing 
processes associated with stellar turbulence, arguably 
the greatest source of uncertainty in progenitor struc- 
ture, besides mass loss, at the time of core collapse. An 
effort to use precision observations of stellar parameters 
to constrain theoretical modeling is also described. 
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1 Introduction 

To first order, stellar evolution theory is a success: an 
evolutionary picture unifying the gross statistical prop- 
erties of stars and stellar populations has been devel- 
oped. However, a number of current astrophysieal puz- 
zles demand a more accurate theory of star evolution 
than currently exists. Examples include (but are in no 
way restricted to): (1) understanding the initial con- 
ditions for explosions in massive stars, including core 
collapse supernovae and gamma ray bursts, (2) the 
formation of 0-Ne-Mg cores and the possibility of ac- 
cretion induced collapse (AIC), (3) the nucleosynthesis 
in AGB stars, including heavy element (Z>26) forma- 
tion through the slow neutron capture process, or s- 
process, and (4) the internal structure of the Sun, and 
the tachocline in particular. All of these research areas 
rely heavily on stellar evolution calculations and suf- 
fer from unaeceptably large theoretical uncertainties, 
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in most cases related to an inadequate treatment of 
turbulent mixing processes. 



2 The current framework for computing 
evolution 

A ID Model. The basic framework for stellar evolu- 
tion relies on the approximation of spherical symmetry. 
When rotation is strong enough to distort the star's 
shape from spherical symmetry the radial coordinate 
is replaced by an equipotential coordinate and quan- 
tities averaged over equipotential surfaces are evolved 
instead, 
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The hydrodynamics governing stellar interiors are 
grossly simplified from the outset by this reduction in 
dimensionality, and a quasi-static approximation is gen- 
erally adopted. Any additional, non-radial hydro- and 
magneto-dynamic processes (such as convection and in- 
ternal waves) are then incorporated on top of this quasi- 
static, ID evolution in the form of transport processes, 
which will be discussed further below. The quasi-static 
equations for single star stru cture and evolution (e.g. 
Kippenhahn fc Weigert 199Clf ) involve the ID formula- 
tions of mass, momentum, and energy conservation to- 
gether with the auxiliary equations 
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which describe the mode of energy transport (eq. [T]) 
and the compositional evolution of the stellar plasma 
(eq. [H, where standard notation has been used and 
the temperature gradient is written in terms of V = 
d\nT/d\nP (not to be confused with the divergence 
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operator in eq. [2]) . The compositional flux is due to 
the cumulative effects of a variety of processes including 
turbulence, internal wave motions, and radiative levi- 
tation and elemental settling, while Xf"'^ is the change 
in composition due to nuclear burning. 

The time evolution of the star is driven by the nu- 
clear transmutation of matter and its heating, energy 
losses at the surface, neutrino cooling in the interior, 
and the indelible mixing of composition and angular 
momentum. The theory of stellar structure and evolu- 
tion consists of trying to capture the dynamics of this 
complex system with fidelity. Additional equations are 
required to incorporate the effects of angular momen- 
tum and magnetic fields, as well as mass loss from the 
surface and binary interactions, including tidal distor- 
tion and mass transfer. 

Energy Transport and Convection. But even when 
considering non- rotating, non- magnetic, single star 
evolution the theory remains incomplete. Energy trans- 
port and thermal instabilities in the stellar interior, and 
their attendant mixing processes, remain some of the 
most pernicious outstanding issues. Energy is carried 
by radiation, thermal convection, or a combination of 
both. The theories underlying these different modes of 
energy transport determine the equilibrium tempera- 
ture gradient V appearing in eq. [T] Radiation transport 
is simplified in the deep interior where photon mean free 
paths (m.f.p.s) are small (compared to any radial scale 
height in the star) and a diffusion approximation is ac- 
curate. In this case the temperature structure depends 
primarily on the opacity of the plasma k and one finds 

V Vrad = 3KLP/{16TTacGmT'^). 

When Vrad > Vad = {d\nT/d\nP)s the stellar 
plasma is unstable to thermal convection and a hy- 
drodynamic flow ensues and modifies the underlying 
temperature stratification. In these regions a model 
for convective energy transport is required to calculate 
the temperature (or entropy) profile of the star. The 
almost universally adopted model is mixing length the- 
ory (MLT) which provides a functional dependence for 

V on the stellar luminosity, structure, and microphys- 
ical properties of the plasma (e.g., opacities, equation 
of state, nuclear burning rates) which is only slightly 
more complicated than the expression for Vrad defined 
above. 

Finally, what is arguably the most complex sit- 
uation regarding energy transport is when the con- 
vective instability extends into the photosphere in 
which case the diffusion approximation breaks down 
and one needs to treat simultaneously thermal convec- 
tion and a complex radiation transport problem. In 
this case as well, the almost universally adopted ap- 
proach is to use MLT and radiative diffusion in evo- 



lution m odels, despite the obvious physical shortcom- 
ings (see Arnett et al.ll201ol for a discussion of the solar 
case). Three dimensional radiation- hydrodynamic sim- 



ulations (e.g.,, Nord lund 1982: Steffen fc Frevtad 11991 
Stein fc Nordlundlll998HRobinson et al.ll2003l ) are be 



ginning to provide self-consistent solutions to this com- 
plex problem, facilitating detailed comparison between 
model and observation. A stro nger foundation for inter 



preting helios eismic data (e.g. IStein fc Nordlundll2001 



Li et al.ll2002f ) and inferring the abunda nces of various 



elemental constituents o f the sun (e.g. lAsplund et al 



2005 : Caffau et al. 20081) are benefitting tremendously 
from this work: these almost brute force approaches to 
solar surface convection are offering sufficiently precise 
information about interior structure to shed light on 
deeper astrophysical con undrums including the f orma- 
tion scenario of the sun ()Guzik fc Mussackll2010l) and, 
in combination with solar neutrino obser vations , stan- 
dard model physics ( Serenelli et al. 20091) . 



Other Instabilities and Mixing. Instabilities involv- 



ferential rotation 



netic fields (jSpruit 



ing composition gradients (ILanger et al.l Il983t ISpruit 
1992t lEggleton et al] boodlMocak et al.l boiObI ). dif- 



insonneault et al. 19891 ) . aiid mag- 
2OO2I) are just as important as ther- 



mal convection with regards to interior mixing. This 
mixing is almost invariably treated as a diffusive pro- 
cess, taking the divergence of the compositional flux to 
be 
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The Eulerian diffusion coefficient D = {dm/dr)~^D is 
most often formulated from a combination of charac- 
teristic length I, velocity v, and tim e scale r, com- 
puted from linear stabih ty theory (jChandrasekhar 
196ll: brazin fc Reidll2004l) . with - or D - P/t. 
Heger et al.l |200d 120051 ) provides a fairly comprehen- 
sive presentation of stellar physics in this "diffusive" 
spirit. Caution is advised when using these results, 
however, since the development and mixing rates for 
these instabilities are not necessarily being represented 
correctly by this approach, particularly in light of the 
fact that the high Reynolds number characterizing 
these flows ensure that the subsequent development 
will be highly turbulen t and in the non-linear regime. 
See also iBruenij (|2005l ) for a related discussion on the 
modern stellar evolution model and its shortcomings in 
core-collapse progenitor modeling. 

Two big open question are: (1) Is a diffusive treat- 
ment for mixing (as in eq. [3]) appropriate? (2) If not, 
how shall we improve this situation? Two processes 
have already been shown to differ dramatically from 
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Fig. 1 Overshoot parameter estimates from wide eclips- 
ing binary (diamonds) and asterose ismological da t a (tri - 
angle s). The values a dopte d by iGirardi et alj (|2000l l 
and IPietrinferni et aD (|2004h for populations synthe- 
sis are labeled G2000 and P2004. The mass un- 
certainties for the eclipsing binary data is negligible 
on the log(M) scale used. Oversh oot parameters in- 



ferred from asteroseismic data are from Soriano fc VauclaiJ 



pOldi). lOi Mauro et id] (|2003|). ISuafez et al.l (120091), 
iBriauet et al.l (l2007f) JPupret et a l.' ("200% lAerts et all 



( 20061 '). iDesmet et al.l (|2009l ') . andMazunida r et al.l 



and t he values from eclipsing binary data are from IClaretl 
(|2007l ). 



diffusion, including the flux of kinetic energy and the 
erosion of convectivc boundary layers due to convec- 
tive overshoot. The kinetic energy flux is now known 
to carry a non-negligible fraction of the total luminos- 
ity of the star under many circumstances, contrary to 
MLT which takes this flux to be zero. In the rare 
case that this contribution to the luminosity is included 
in a stellar model it is usually treated using a "down 
gradi ent approximation" with Fj^ oc -VEk iKuhfusd 
19861) ■ which is in direct conflict with simulated flow 



as well as basic conside rations of energy conservation 
( Meakin fc ArnettI boiol) . Similarly, the evolution of 
convective boundaries are likely to evolve in a manner 
which is far from diffusive in nature, and behave instead 
as sharp entrainment interfaces, which will be discussed 
further in [J5l 



3 How well does the current theory work? 

Quantifying Uncertainty. While the uncertainty inher- 
ent in the stellar modeling framework is difficult to 
quantify, a formal lower bound can be assessed by com- 
paring stellar parameters calculated in the following 
ways. (1) The results from a variety of evolution codes 



Fig. 2 The ages and masses of wide eclipsing binary 
star data having p recision observational parameters from 
iTorres et all (j2010[ ). These are all upper main sequence 
stars with convective cores and radiative envelopes and are 
therefore some of the simplest stellar objects to model. The 
dashed curve represents an approximate main sequence turn 
off age as a function of initial mass. 



which have nominally the same input physics and pa- 
rameters can be compared. This exercise addresses the 
question: How accurately is the mathematical model 
being captured by the numerical implementation? (2) 
The results from a single code using a range of plau- 
sible adjustable parameters and input physics can be 
compared. This addresses the question: How sensitive 
is the stellar model to uncertainties in the input physics 
and parameters within the framework being used? 

While some effort has been exerted to obtain con- 
sistency between di fferent implementations of the same 
physi cal model (e.g. lMontalban et al ][2nn8HMarconi et ahj 



20081 ). code to code scatter persists even for the rela- 



tively simple case of main sequence evolution. A recent 
comparison between the inferre d properties of ma ssive 



O and B stars within 500 pc (IHohle et al.l 120101) us 



ing three modern evolution code s (jSchaller et al 



1992 



Bertelli et al.lll994t IClaretll2004l ) indicate uncertainties 
at the '^20 - 25 % level. Similar degrees of uncer- 
tainty were found when computing the mass limits as- 
sociated wi th 0-Ne-Mg core forma tion in intermediate 



mass stars (jPoelarends et al.1120081 ). 



Of deeper concern is the range of evolutionary out- 
comes found for the same initial conditions using what 
are apparently an acceptable range of plausible input 
physics and adjustable parameters. This is discussed 
next. 

Mass Limits. An informative way to summarize a 
stellar evolution model is to calculate various mass lim- 
its. Examples include: The mass above which a star will 
form a black hole, Mbh; the mass above which core he- 
lium burning ignites under non-degenerate conditions. 
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Mfjo-f; the mass above which an iron core forms and 
undergoes gravitational cohapse, Mpc-cc; etc. A sum- 
mary of such Hmits for massive star evolution, with 
an empha sis on corc-c o Uapse supernova theory, is pre- 
sented by iHeger et alj (|2003l ) ; and the fate of massive 
AGB stars and the potential fo r -Ne-Mg core col- 



lapse a re discussed bv lSiessI (|2007l ) and lPoelarends et al 
( 20081) . 



By comparing the mass limits found from a finely 
enough sampled grid of models and observational stel- 
lar data, the theoretical model can be tested. The 



landmark paper by iMaeder fc MermilliodI (|l98lh used 
this approach to examine the consistency between stel- 
lar cluster data and their evolution model, and con- 
cluded that additional mixing can bring observation- 
ally inferred and theoretically calculated mass limits 
into better agreement, and summarized this in terms 
of an "overshoot" parameter ~ 0.2, a quantitative 
value which remains preferred among stellar population 
modelers f or stars having masses larger than a few time s 
solar (e.g 



3r stars having masses larger than a tew time s 
Girardi et al.ll2000l : IPietrinferni et al.ll2004l ). 



In addition to the overshoot parameter, which is a 
measure of how much material mixes beyond the lim- 
its of a formally defined unstable region, the mixing 
processes mentioned in ^ also shift mass limits. The 
double diffusive instability arising in thermally unstable 
regions which are stabilize d by compositio n gradients, 
dubbed " semiconvection" (jMerrvfieldll 1 9951 ) . is particu- 
larly important in massive star evolution. This process 
can change the mass of the He core following main se- 
quence evolution, perhaps the parameter most indica- 
tive of how a massi ve star will end its life, by a factor 
of ^ 50% or more (jWooslev fc Weaver! Il988l : IStaritsin 
20091) . When to include this mixing process or not, and 
its strength s when included, remain open questions (e.g. 
BiellolboOll) . 

Fortunately, interesting constraints can be placed 
on stellar evolution theory by observational data when 
the observational uncertainties are less than a few per- 
cent. Both asteroseismic and wide eclipsing binary 
data, which we discuss next, are beginning to meet 
these requirements. 



4 The need for precision observational data 

Wide Eclipsing Binaries. Observationally determined 
effective temperatures (Ti,T2), luminosities (Li,L2), 
radii (i?i,i?2), and masses (Mi,M2) are known for 
many of these sys tems to better than a few percent 
( Torres et al1l2010f l. In some cases precision surface ro- 
tational velocities ([i;sini]i, [wsini]2) and compositions 
Xi^i — Xi^2 are also available. A stellar evolution model 




Fig. 3 A detailed view of the boundary layer separating 
a turbulent convection zone from an over lying: stably strat- 
ificd layer in a stellar interior simulation (|Meakin fc ArnettI 
[2007b, 20lJ). The stellar radius increases vertically upwards 
in the right-most panel while the other four panels are hor- 
izontal slices having radial positions indicated by the black 
lines in the right-most panel. The grey scale represents oxy- 
gen abundance, with lighter values indicating higher abun- 
dance. 



can be tested against this data by comparing, e.g., 
the model radii and temperatures i?2; Ti, T2) = 
/(Ml, M2, Xi] iage, Q!j), where the hat symbol indicates 
model data. Here the stellar masses and compositions 
arc also known observationally to high precision. The 
age of the binary system tagc is a fitting parameter, 
and the theoretical model is repre sented by th e func - 



Ribasetal 



|200fll ) 



tion / and the p arameter set Uj. 
and IClareti (|2007l ) have studied the dependence of the 
convective overshoot parameter ctov on the stellar mass 
using this type of procedure. 

A stero seismology. The set of normal mode oscilla- 
tions frequencies {vk} found by time monitoring stel- 
lar luminosity, together with spectroscopic (log 5 and 
Toff) and photometric information (e.g., parallax tt 
and luminosity L), provides another important test 
of stellar evolution theory. In this case, the model 
data {i>fc, Tcff , log^, L} = /(M, X^, iago; aj) is compared 
to the observed data to find a best fit stellar mass, 
compositioj^, age, a nd model parameters aj (see e.g. 



Vauclair et al.l 120081 ). The strength of this method lies 
in the large number of observable frequencies, each hav- 
ing a unique spat ial dependence on internal structure 



( Unno et al.lll989l) 



Results and Implications. A summary of the over- 
shoot parameter and its dependence on stellar mass, 
inferred from wide eclipsing binary and asteroseismo- 
logical data, is presented in Fig. [TJ While the error 
bars are still quite large in Fig. [U the scatter in the 
data may be indicating that we are missing some essen- 
tial physics. A simple and intuitive possibility is that 
"overshoot" might not just have a mass dependence but 
also a time dependence. We are just beginning to col- 
lect enough precision data to investigate questions like 



^Observed values of [Fe/H] can be used to estimate the helium 
abundance Y by adopting a cosmic enrichment rate. Uncertain- 
ties and intrinsic scatter in these laws generally undermine con- 
fidence in this approach so that Y is generally left loft as a free 
parameter to be fit. 
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this. In Figure [2] we present a subset of wide eclipsing 
binary data for upper main sequence stars in the mass-^ 
age plane, showing that we can now study stars of a 
given mass at multiple points in time during the course 
of their evolution (Meakin, Sukhbold, & Arnett, 2010 
in preparation). 



5 Fiiture directions 



can be used to test mixing laws like the one described 
here, which has a natural time dependence and physical 
consistency built into it. This type of theoretical mod- 
cling provides an essential bridge between turbulence 
simulation data, like that shown in Fig. [5] whic h is used 
to dev elop these models, and reality . (See also I Meakin 
( 20081) and lArnett fc MeakinI (|201ol ) for additional dis- 



cussion concerning the use of simulation data to study 
stellar turbulence.) 



Precision observational data is beginning to reveal 
deficiencies in our theoretical models, and improve- 
ments are needed. More highly parameterized mod- 
els offer one means to accommodate new data, but 
rarely lack predictive power and even less frequently 
motivate new discoveries. A more satisfying ap- 
proach is to capitalize on what the new data can 
tell us by pursuing a more comprehensive theoreti- 
cal picture. Fortunately, many promising develop- 
ments in understanding turbulent flow have already 
taken place in related fields of study w hich await ap- 
plication in stellar evolution (see e.g. ICanutd Il992 . 
and references therein) in addition to recent work 
dedicated specifical ly to understand i ng stellar inte- 



rior turbulence (e.g. Bazan et al 



2004: Rogers fc Glatzmaier 2006: Eggleton et al 



Herwig et al.ll2006 



200i iBrowning et al 



20M,|2007b; 



2006 



6 Conclusions 

Stellar evolution, particularly our understanding of 
massive stars, is a crucial input to many areas of astro- 
physical research. Our current best models, however, 
suffer severe deficiencies in treating turbulent trans- 
port and mixing to the degree that it is the domi- 
nant uncertainty (in combination with mass loss) in our 
calculations. These shortcomings are currently being 
addressed by studying ever more realistic simulations 
of turbulent flow and incorporating the results from 
knowledge garnered in related fields of inquiry (e.g., 
geophysical). This work is beginning to provide deep 
insights into stellar interior mixing processes. While 
this work is still in a relatively formative stage, early 



Mocak etrffldi^ s are beginning to reveal promising alternative 



. .Meakin fc Arnettll 
20081 12OO9I l2010aibl : ICaraud et al.ll20l'dir 

One such striking development is the recognition 

that a turbulent convective bounda ry is like ly to 

evolve as an entrainment interface (e.g. Fernandol 1991 ; 

Meakin fc Arnetti2007bl :l Mocak et al.ll2 009l: IWoodward et aic| | ning stellar evolution theory. Fortunately, precision 

data is beginning to arrive just in time to begin test- 
ing the increasingly more sophisticated models of stellar 
interior physics that are presently being formulated. 



approaches to the back of the envelope, mixing-length 
style treatments presently being used. 

Finally, it can not be emphasized enough how im- 
portant precision observational data is to testing and 



20091 ). While this svnergv between geophysics and as- 
trophysics is very intuitive and obvious in retrospect, 
it nevertheless provides a more robust physical basis 
for describing "overshooting" and its time dependence, 
which has utility. From relatively simple physical argu- 
ments once can deduce that the boundary of a turbulent 
convection zone in the deep inter ior should evolve at 
a speed given by w = (T„Ai?ig" (jKantha et al.lll977l ) 
where cTu is the r.m.s. turbulence velocity near the 
boundary, Riq is the bulk Richardson numbeiH, A is an 
" efficiency" parameter which depends on the turbulence 
in the vicinity of the interface, and n is an exponent 
close to 1. The data for the stars included in Fig. [2] 



^The bulk Richardson number is a dimensionless measure of the 
potential energy in stratification separating the turbulent layer 
from the stably stratified layer, often written in terms of the buoy- 
ancy jump Sb and integral length scale of turbulence L, compared 
to the strength of the turbulence (t„, with Rig = 5bL/a^. Small 
values indicate boundaries which are strongly distorted by the 
adjacent turbulence, while large values indicate relatively undis- 
turbed boundary layers. 
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